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We study the charge noise Sq in Josephson qubits produced by fluctuating two level systems 
(TLS) with electric dipole moments in the substrate. The TLS are driven by an alternating electric 
field of angular frequency Q and electric field intensity /. It is not widely appreciated that TLS in 
small qubits can easily be strongly saturated if / 2> I c , where I c is the critical electric field intensity. 
To investigate the effect of saturation on the charge noise, we express the noise spectral density in 
terms of density matrix elements. To determine the dependence of the density matrix elements on 
the ratio I/I c , we find the steady state solution for the density matrix using the Bloch-Redfield 
differential equations. We then obtain an expression for the spectral density of charge fluctuations 
as a function of frequency / and the ratio I/I c - We find 1// charge noise at low frequencies, and 
that the charge noise is white (constant) at high frequencies. Using a flat density of states, we find 
that TLS saturation has no effect on the charge noise at either high or low frequencies. 

PACS numbers: 74.40.+k, 03.65.Yz, 03.67.-a, 85.25.-j 



I. INTRODUCTION 

The superconducting Josephson junction qubit is a 
leading candidate in the design of a quantum computer, 
with several experiments demonstrating single qubit 
preparation, manipulation, and measurement as 
well as the coupling of qubits.— £ A significant advantage 
of this approach is scalability, as these qubits may be 
readily fabricated in large numbers using integrated cir- 
cuit technology. However, noise and decoherence are ma- 
jor obstacles to using superconducting Josephson junc- 
tion qubits to construct quantum computers. Recent 
experiments 7 ^ indicate that a dominant source of deco- 
herence is two level systems (TLS) that are fluctuating 
in the insulating barrier of the tunnel junction as well 
as in the dielectric material used to fabricate the circuit, 
e.g., the substrate. The two level fiuctuators that have 
electric dipole moments can induce image charges in the 
nearby superconductor and hence produce charge noise 
Sq (w) . 9 - 10 - 11 - 12 - 13 - 14 - 15 - 16 

Previous theories of charge nois o 14 ' 15 i 16 have neglected 
the important issue of the saturation of the two level 
systems by electric fields used to manipulate the qubits. 
Dielectric (ultrasonic) experiments on insulating glasses 
at low temperatures have found that when the elec- 
tric (acoustic) field intensity / used to make the mea- 
surements exceeds the critical intensity I c , the dielec- 
tric (ultrasonic) power absorption by the TLS is satu- 
rated, and the attenuation decreases as the field inten- 
sity increases i & 17 ' 18 i 19 i 20 (If £ cos(fM) denotes the electric 
field, then we define the intensity / = £ 2 .) Previous the- 
ories of charge noise in Josephson junctions assumed that 
the TLS were not saturated, i.e., that / -C I c . This seems 
sensible since charge noise experimentsSA have been done 
in the limit where the qubit absorbed only one photon. 

However, the following simple estimate shows that 
stray electric fields associated with this photon could sat- 
urate two level systems in the dielectric substrate which 



supports the qubit. We can estimate the voltage V across 
the capacitor associated with the substrate and ground 
plane beneath a Cooper pair box (see Fig. [T]) by setting 
CV 2 /2 — hu> where huj is the energy of the microwave 
photon. We estimate the capacitance C = e e r A/L ~ 7 
aF using the area A = 40 x 800 nm 2 of the Cooper pair 
box, the thickness L = 400 nm of the substrate^ and the 
relative permittivity e r — 10. Using oj/2tt = f = 10 GHz, 
we obtain a voltage of V ~ 1.4 mV. A substrate thick- 
ness L of 400 nm yields an electric field of £ ~ 3.4 x 10 3 
V/m. 

We can compare this with the critical intensity I c and 
the associated critical electric field £ c which has been 
measured experimentally^. Martinis et al& measured 
the low temperature dielectric loss tangent of amorphous 
Si0 2 at / = 7.2 GHz and amorphous SiN^ at / = 4.7 
GHz as a function of the root-mean-square (rms) volt- 
age. They found that the loss tangent was constant at 
low power, but rolled over and decreased above a critical 
rms voltage V c ~ 0.2 fiV. For a capacitor thickness of 
300 nm, the associated critical field is £ c ~ 0.7 V/m. So 
£/£ c - 5 x 10 3 , and I/I c = (£/£ c f - 2 x 10 7 > 1. 

We can do a similar estimate to show that a single 
photon would even more strongly saturate resonant TLS 
in the insulating barrier of the tunnel junction. We use 
the same numbers as before but with C = 1 £F and the 
thickness of the junction L — 1.5 nm to obtain £ ~ 7xl0 4 
V, and I/I c ~ 10 10 > 1. 

However, there are only a few TLS in the oxide bar- 
rier of a small tunnel junction as the following estimate 
shows. For a parallel plate capacitor with L = 1.5 nm and 
A = 1 ^im 2 . the volume is V Q = 1.5 x 10 -21 m 3 . The typ- 
ical TLS density of states^ 2 - is 10 46 /Jm 3 . However, only 
a fraction of these have electric dipole moments. So we 
will assume that the density of states of TLS with electric 
dipole moments is Ptls — 10 45 / Jm 3 ~ 663//iGHz/im 3 . 
Using this value of Ptls, we find that in a small tun- 
nel junction, there are only 2 TLS with an energy split- 
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FIG. 1: (Color online) Distribution of two level systems (TLS) 
in the dielectric substrate of a Cooper pair box. The super- 
conducting island with the Josephson junction qubit is the 
black rectangle on the top of the substrate. There is a ground 
plane underneath the substrate. 

ting less than 10 GHz. A single fluctuates would have 
a Lorentzian noise spectrum. The presence of 1// noise 
implies many more than 2 fluctuators. It is likely that 
these additional fluctuators are in the substrate. Our 
main point is that TLS in small devices are easily satu- 
rated. It is therefore important to analyze the effect of 
TLS saturation on the charge noise both at low and high 
frequencies / of the noise spectrum. 

In this article, we explore the consequences of this sat- 
uration on the spectral density of polarization and charge 
fluctuations. We consider a driven system consisting of 
two level systems with electric dipole moments that fluc- 
tuate randomly leading to fluctuations SP(t) in the po- 
larization denoted by P(t). In addition, the dipole mo- 
ments of these TLS couple to an applied ac electric field 
that drives the system with an angular frequency f2. 

Let us consider a single two level system. According to 
the Wiener-Khintchine theorem, in the stationary state 
the polarization noise spectral density Sp(lu) of this two 
level system is twice the Fourier transform of the polar- 
ization autocorrelation function: 

/+oo 
d(t 1 -t 2 )e^- t ^[{SP(t 1 )SP(t 2 ))} p 
-oo 

/+oo 
d(ti - hy^-^TripiSPit^SPih))} 
-oo 

where [...] p denotes a trace over the density matrix p, and 
(...) denotes an average over the time series. 5P(t) = 
[P(t) — (P)} is the fluctuation in the polarization P(t) at 
time t and (P) is the polarization averaged over the time 
series. From the convolution theorem, Sp(oj) ~ |5P(w)| 2 . 
If the average over the density matrix is equivalent to an 
average over the time series, then one can easily show 
that SP(u) = 0) = and hence, S p (uj = 0) = 0. 

We assume that the total density matrix px is the 
product of two factors: one factor p that contains the 
contributions of the ac driving field and the other that 
contains the contributions of the random fluctuations 
of the dipoles. We make this division because we can 



find the time dependence of the density matrix (in the 
Schroedinger representation) due to the driving field by 
solving the Bloch-Redfield equations. Since it is not clear 
how to include the random fluctuations of the two level 
systems in the Hamiltonian and hence, in the Bloch- 
Redfield equations, we will treat the random fluctuations 
separately So the density matrix p contains the contri- 
butions of the driving field to Sp(uj) while the polariza- 
tion autocorrelation function (5P(ti)5P(t2)) contains the 
contributions of the random fluctuations of the electric 
dipoles. We will solve the Bloch-Redfield equations in 
steady state to find the time evolution of the density ma- 
trix p{t) and its dependence on the ratio I/I c (Cl, T) of the 
electric field intensity / to the critical intensity 7 c (f2,T) 
which is a function of the driving angular frequency fl 
and the temperature T. The benefit of this approach 
is that it is valid in both equilibrium as well as steady 
state non-equilibrium situations. From the polarization 
noise spectral density, we can obtain the charge noise 
Sq(u). We then average over the distribution of inde- 
pendent TLS. Unlike previous theoretical efforts , 14 i 15 i 16 
we use the standard TLS density of states that is a con- 
stant independent of energy. 

At low frequencies [hf -C kpT) the system is in equi- 
librium, and we find 1// charge noise that is propor- 
tional to the temperature and to the dielectric loss tan- 
gent tan 5 that has a well known contribution from the 
electric dipole moments of TLS j 17 ' 22 ' 23 In addition the 
low frequency charge noise has negligible dependence on 
the electric field intensity ratio I/I c - 

At high frequencies (hf ^> kpT) we find that the 
charge noise is white noise independent of frequency. It 
has a very weak dependence on the ratio I /I C (£L, T) and 
the driving frequency fi. We also find that the amplitude 
of the high frequency white charge noise decreases grad- 
ually as the temperature increases. The fact that the 
charge noise spectrum depends very weakly on the ra- 
tio i"/J c (Q, T) indicates that the saturation of two level 
systems does not affect charge noise. 

The paper is organized as follows. In Section [Til we 
present our model of a TLS in an external driving field. 
In Section IIII1 we use the fluctuation-dissipation theo- 
rem to give an expression for the charge noise in ther- 
mal equilibrium. In Section IIV( we take a more general 
approach that is valid in both equilibrium and nonequi- 
librium cases. In particular, we derive a general ana- 
lytic expression for the spectral density of polarization 
and charge fluctuations of an individual two level system 
(also referred to as a fluctuator) in terms of the density 
matrix. In Section [V] we solve the Bloch-Redfield lin- 
ear differential equations for the density matrix. We find 
the steady state solution of the Bloch-Redfield equations 
and we analyze its dependence on the ratio I / I C (VL,T). 
In Section IVTl we investigate the noise spectrum of a sin- 
gle random telegraph fluctuator. We then average over 
the distribution of independent TLS numerically to de- 
termine the frequency dependence of the noise spectrum. 
A summary is given in Section [VIII 
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II. TWO LEVEL SYSTEM (TLS) 

In applying the standard model of two level systems 
to Josephson junction devices, we consider a TLS that 
sits in the insulating substrate or in the tunnel barrier, 
and has an electric dipole moment p consisting of a pair 
of opposite charges separated by a distance d. The elec- 
trodes positioned at z = and z — L are kept at the 
same potential. The angle between p and the z-axis, 
perpendicular to the electrodes, is 9. The dipole induces 
charge on the electrodes. As we show in the appendix, 
the magnitude of the induced charge Q on each electrode 
is proportional to the z-component of the dipole moment, 
p z = pcos9, i.e., 



pcost 



(1) 



The dipole flips and induces polarization fluctuations and 
hence charge fluctuations on the electrodes. 

The TLS is in a double-well potential with a tunneling 
matrix element Ao and an asymmetry energy A.— The 
Hamiltonian of a TLS in an external ac field can be writ- 
ten as H(t) = H + Hi(t), where H = \{Ao z + A <r x ), 
and Hi(t) = —cr z p ■ £ a c(t)- Here a x ,z are the Pauli spin 
matrices and £ oc (i) = £ ac cos Qt is a small perturbing ac 
electric field of angular frequency that couples to the 
TLS electric dipole moment. £ oc points along the z— axis. 

After diagonalization of Hq, the Hamiltonian becomes 



H(t)=H + H 1 {t) 

Ho = \ev z 
H\{t) — —rj(A a x + Acr z )cosfit, 



(2) 

(3) 
(4) 



where E = yj A 2 + A 2 D is the TLS energy splitting and 
V = P ' £ac/E- Notice that 77 is a small dimensionless 
variable (77 w 5 x 1CT 3 for p = 3.7 D, £ ac w 10 3 V/m 
and E/h * 10 GHz). (The dipole moment of an OH - 
molecule is 3.7 DM) The energy eigenbasis is denoted by 
{|+), |— )}, and the corresponding eigenvalues are E± = 
±E/2, where + (— ) refers to the upper (lower) level of 
the TLS. The energy splitting E will also be referred to 
as hujQ. 

An excited two level system can decay to the ground 
state by emitting a phonon. The longitudinal relaxation 
rate is given by>^ 



EAi coth 



E 



2k B T 



(5) 

where p is the mass density, eg is the longitudinal speed 
of sound, Ct is the transverse speed of sound, and 7^ is the 
deformation potential. In this paper we will use the val- 
ues for Si0 2 : 7d =1 eV, p = 2200 kg/m 3 , q=5800 m/s, 
and c t =3800 m/s. Typically, T\ varies between 10 -9 s 
and 10 4 s for temperatures around 0.1 K. The distribu- 



tion of TLS parameters is given b y 22 > 25 

Ptls 



P(E,T 1 ) = 



2T l ^l-{T mm {E)/T l ) 



(0) 



where Ptls is a constant density of states that represents 
the number of TLS per unit energy and unit volume. The 
minimum relaxation time T m in (E) corresponds to T\ for 
a symmetric double-well potential (i.e., E — Aq): 



ll 



(2irph 4 ) 



1 



E A coth 



E 



2k B T 



(7) 

Alternatively, the TLS distribution function can be ex- 
pressed in terms of the TLS matrix elements A and Aq: 



P(A,A 



Ptls 
Aq ' 



(8) 



The typical range of values for A and Ao are < A/fee < 
4 K and 2 pK < A /k B < 4 K, where k B is the Boltz- 
mann's constant. We will use these values for our numer- 
ical integrations in Section |VT1 



III. THERMAL EQUILIBRIUM EXPRESSION 
FOR CHARGE NOISE 



We begin by considering the case of thermal equilib- 
rium. According to the Wiener-Khintchine theorem, the 
charge spectral density Sq(u>) is twice the Fourier trans- 
form ^q(uj) of the autocorrelation function of the fluc- 
tuations in the charge. In equilibrium we can use the 
fluctuation-dissipation theorem 2 ^ to find that the (un- 
symmetrized) charge noise is given by: 



S Q (k,Lj) = 



4ft 



^ Q — htu/k£ 



(9) 



where Q is the induced (bound) charge on the electrodes 
and Xg(fc,a;) is the Fourier transform of 

x5(r.^,O- " 0(r - t) -g''- OU - do) 

where [...] c is a commutator, and (...) e is an ensemble 
average. We use Q = J P ■ dA, where P is the electric 
polarization density, and choose P z and gL4||z since Q ~ 
\p z \ to find 



X '! i {k,u)=e A 2 ^{k,u) 



ill) 



where Eq is the vacuum permittivity, A is the area of 
a plate of the parallel plate capacitor with capacitance 
C, and x'p s , oS) is the imaginary part of the electric 
susceptibility. We set k = 0, and use 



£ Xp z (i^) = e'(uj) tan 5 (w) 



(12) 



where the dielectric loss tangent tan(5(a>) = e" (w) / e' (u>) . 
s'(w) and e"(w) are the real and imaginary parts of the 
dielectric permittivity, respectively. We also use 



C = e'A/L 



(13) 
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to find 

where Sq(lu) = Sq(k = 0,u>)/V o , the volume of the ca- 
pacitor is V = AL, and e'(ui) = e' + etls{ uj ) — e' = Eo£ r 
where e r is the relative permittivity. The frequency de- 
pendent permittivity £tls(w) produced by TLS is neg- 
ligible compared to the constant permittivity e' The 
TLS dynamic electric susceptibilities (x'(w), x"( w ))7 an d 
hence the dielectric loss tangent, can be obtained by solv- 
ing the Bloch equations in cquilibriumj 17 i 19 i 27 One can 
then average over the distribution of TLS parameters. 
However, since we will be considering driven systems that 
are in a noncquilibrium steady state, we need to take a 
more general approach which is described in the next 
section. 



IV. GENERAL EXPRESSION FOR SPECTRAL 
DENSITY OF POLARIZATION AND CHARGE 
FLUCTUATIONS OF A TWO LEVEL SYSTEM 

The noise in our model is due to a fluctuating two level 
system with an electric dipole moment that changes its 
orientation with respect to the direction of the applied 
driving field while keeping its magnitude constant. In 
this section, we begin by deriving a general expression for 
the polarization noise Sp(co) of a single TLS that is valid 
at all frequencies and in both equilibrium and nonequi- 
librium situations. Since we are interested in TLS satu- 
ration, this formulation will apply to a driven system in 
nonequilibrium steady state. We then relate the polar- 
ization noise to the charge noise Sq(lv). 

According to the Wiener-Khintchinc theorem, the po- 
larization spectral density Sp(u>) in the stationary state 
is twice the Fourier transform of the autocorrelation func- 
tion of the fluctuations in the polarization: 



/+oo 
d(ii - t 2 )e lu{tl - t2) 
-oo 

x Tt{ph [(Pff(ti)P ff (t 2 )> - (P) 2 ] }, (16) 

where the density matrix pn is time independent in the 
Heisenberg representation. In the Schrodinger represen- 
tation, the density matrix has time dependence. Wc 
now change from the Heisenberg representation to the 
Schrodinger representation (denoted by the subscript S) . 
Recall that an operator A nit) in the Heisenberg repre- 
sentation can be expressed in the Schrodinger represen- 
tation by A H (t) = W(t,t )As(t)U(t,t ), where U(t,t ) 



is the unitary time evolution operator. Hence 

Tr{p H PH(ti)P H (t2)} = TKAs(ti)Ps(t0l7(ti,t2) 

= F(t u t 2 ). (17) 

To simplify the notation we have temporarily omitted 
the symbol (...) denoting the time average. The spectral 
density of polarization fluctuations in the Schrodinger 
representation becomes: 

/+oo 
d(ti - t^e^-^FiU, t 2 ) - 2{P) 2 S(lu), 
-oo 

(18) 

where we assume that the system is in a stationary state 
so that the function F(t\,t2) depends on (ti — fa). It can 
be expressed as 

F{t u t 2 )= <"»(*i)|p(*i)|n(*i))M*i)l^(*i)lP(*i)> 

m,n,p 

x (p(ti)|t7(t x ,t 2 )P(t2)l7(t2,ti)|m(ti)) 

- ^ (mihMtxMhMnih^PihMh)) 

m,n,p 

x (p(t 2 )\P(t 2 )\m(t 2 )), (19) 

where p(t) is the density matrix in the Schrodinger rep- 
resentation, and m, n, and p denote eigenstates of Ho- 

As we mentioned in the introduction, we are consid- 
ering the density matrix p of a single TLS that contains 
the time dependence of the external driving field. The 
random dipole fluctuations are contained in Pjj(t). Let 
a stand for m, n, or p. Then \a(t)) = exp(—iE a t/H)\a) 
and H \a) = E a \a). We now switch from \n(t)}, \m(t)), 
and \p(t)) to the |+) and — ) eigenstates of a TLS to 
obtain: 

F{t u t 2 ) = Wi)P{ti)] + +P++{t 2 ) 
+ [ /9 (t 1 )P(i 1 )]__P__(i 2 ) 

+ e+ <w °( il - t2 )[p(*i)P(*i)]+-P-+(*2), (20) 

ire P aa > (t) denotes the aa'th element of the P(t) ma- 
, [p(t)P(t)) aa > represents the aa'th element of the 
P(t) matrix, and luq = E/h. We will see in Section IVl 

that p-\ and p |_ arc first order in the small parameter 

X] = p • £ac/E. -q <C 1 for both small and large values of 

Ijl c . So we will neglect terms with p_| and p |_. This 

leads to an approximate expression for F{ti,t 2 ): 

F(t u t 2 ) « /9++ P ++ (i 1 )P ++ (i 2 ) 
+ P— P— (ti)P— (t 2 ) 
+ e- <w °(*i-*»)p__P_ + (t 1 )P + _(ia) 
+ e+^-^p ++ P + _(t 1 )P^+(h). (21) 

Let P\\(t) be the projection along the ac external field 
of the polarization operator associated with the dipole 



/+oo 
dfa-tje^-V [(PH(h)PH(t 2 )) - (P) 2 ] p , 
-OO 

where the subscript H denotes the Heisenberg represen- 
tation and (P) is the time averaged value which is in- y 
dependent of the actual representation. We can rewrite t 
this expression in the Heisenberg representation as c 
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moment p of a two level system. P\\(i) has stochastic 
fluctuations due to the fact that the electric dipole mo- 
ment of the two level system randomly changes its orien- 
tation angle 9(t) with respect to the applied electric field. 
Hence, in the TLS energy eigenbasis, we can write 



pcos(fl(t)) fA Q A 



(22) 



where Po(t) = —pcos(9(t))/V and V is volume. 

Substituting P|| for P in Eq. (0 and using Eq. ([2"Tj). 
we obtain 



SpJuj) 



d(ii-t 2 )e^ tl -* 2 )(P (t 1 )P (t 2 )) 



Ao 

E 



e iwo(*i-*2)p_ 



ft-o(ti-ta) p++ ]l _2(P) 2 %), 



(23) 



Since for stationary processes the correlator 
(Po(ii)Po(^2)} is a function of (ii — £2), we can de- 
fine S Po (lu) = 2/+^d(t 1 - ^e^"^ (P (t 1 )Po(i 2 )). 
Then we have 

^P||(w) = [^j S Po( UJ ) 

+ (^) ^ — Sp °( u ~ ^ + P++Sp °^ u> + 

-2{P)H{u J ). (24) 

This is a general formula for the spectral density of the 
polarization fluctuations assuming that the fluctuations 
in the orientations of the electric dipole moments of TLS 
are a stationary process. The last term ensures that 
Sp,, (lu = 0) = 0. The first term, which is proportional to 
Sp (lj), is the relaxation (REL) contribution. It is as- 
sociated with the TLS pseudospin a z whose expectation 
value is proportional to the population difference between 
the two levels of the TLS. The relaxation contribution to 
phonon or photon attenuation is due to the modulation 
of the TLS energy splitting E by the incident photons 
which have energy hio <C E. This modulation causes the 
population of the TLS energy levels to readjust which 
consumes energy and leads to attenuation of the incident 
electromagnetic flux. Because p++ + p = 1, the re- 
laxation term has no dependence on the density matrix, 
and so will not be affected by saturation effects. Since, 
as we will see in section IvTI this term dominates at low 
frequencies, this implies that the low frequency noise will 
not be affected by TLS saturation. 

The middle two terms in Eq. ([24]) are proportional to 
Sp {ujiu!o) and are resonance (RES) contributions. The 
resonance terms are associated with the x and y com- 
ponents of the TLS pseudospin that describe transitions 



between energy levels. They describe the resonant ab- 
sorption by TLS of photons or phonons with hu> = E. 
We will see in section IVII that the resonance contribu- 
tions are dominant at high frequencies. 

To obtain the charge noise Sq(lo) from the polariza- 
tion noise, we make use of the following formulas. In a 
polarized medium, the induced (bound) charge is 



Q 



P-dA 



where P is the electric polarization. We choose P z 
and dA||z since Q ~ \p z \ = \pcos8\. Then 

Sq(u) = A 2 SpXlo) 



(25) 
= P| 

(26) 
(27) 



where A is the area of a plate of a parallel plate capacitor 
with capacitance C = e'A/L, and e' is the real part of 
the dielectric permittivity. 

In this section we have derived expressions for the po- 
larization and charge noise of a single TLS in terms of 
the density matrix. In the next section we will solve the 
Bloch-Redfield equations for the time dependent density 
matrix of a TLS subjected to an external ac driving field. 



V. THE BLOCH-REDFIELD EQUATIONS 



From Eqs. (|T8|) and (jl9|) , we see that we need the time 
dependent density matrix to calculate the polarization 
noise spectrum. In this section, we solve the Bloch- 
Redfield equations to find the time evolution of the den- 
sity matrix of a single TLS subject to an external ac 
electric field. These equations combine the equation of 
motion of the density matrix with time-dependent per- 
turbation theory, taking into account the relaxation and 
dephasing of TLS. 

We follow Slichter— and write the following set of lin- 
ear differential equations for the density matrix elements 

Paa'(t): 



dp 



dt 



P,/3' 



(28) 



where a and (3 can be either + or — , corresponding to 
the energy eigenstates of the TLS, and Raa'.pp' are the 
Bloch-Redfield tensor components which are constant in 
time. They are related to the longitudinal and transverse 
relaxation times, T\ and T2. In Eq. (f28|) . Ho and H\{t) 
are Hamiltonians given by Eqs. © and (g]), respectively. 
The thermal equilibrium value of the density matrix is 
denoted by p eq {T). 

In thermal equilibrium only the diagonal elements 
of the density matrix are nonzero, and are given 
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by pt q _(T) = exp(+E/2k B T)/Z and p%+(T) = 
exp(— E/2kgT)/Z , where the partition function Z = 
[exp(-E/2k B T) + exp(+E/2k B T)}. Here we are using 
the fact that in thermal equilibrium, the density matrix 
can be represented as 



Peq(T) 



1 



Z{T) 



,-H /k B T 



(29) 



One may wonder whether instead we should use 



p(T,t) 



Z{T,t) 



exp 



(.ffo +&(*)) /*bT] (30) 



since the ac field changes the TLS energy splitting. 
Slichter has discussed this issue in his booki 28 Eq. ([29]) 
is appropriate if the TLS are too slow to respond to the 
external ac field, but Eq. ([30| should be used if the exter- 
nal field varies much more slowly than the response time 
of the TLS. In the latter case, the external field looks 
like a static field to the TLS. For the cases of interest, 
typical experimental ac external fields operate at several 
GHz while the response or relaxation time of TLS is T± 
which, as we said earlier, typically varies between 10 -9 s 
and 10 4 s. So it is reasonable to use Eq. (|29j) . 

Eq. ([28]) can be written in the form: 



dp, 



-t{E a i - E a )p aa i + - y"^[/w {a"\H\ (t) \a'} 

a" 

- {a\H x {t)\a")pa»a>] 

+ ^ Raa ' M'\PW - Pp q (3'( T )}- 



(31) 



P,/3' 



Next we use the fact that in the relaxation terms 
Raa', /3/3 1 , the only important terms correspond to 28 i 29 
a — a' — f3 — (3' . In addition, the Bloch-Redfield ten- 
sor is symmetric, so we have the following relations for 
the dominant components: 



R — ,++ 
R-+-+ 



Ra 



l 



R 



+-,+- 



1 



(32) 
(33) 



The longitudinal relaxation time T± is given by Eq. ([5]). 
For the transverse relaxation time T 2 , we will use the 
experimental valued 



8 x 1Q- 7 
T 2 ~ j, sec, 



(34) 



where T is in Kelvin. From relations (|3"2")l and (f3"3")l . the 



set of linear differential equations (|3Tj) becomes 
J-l 



dp 



a/ =^(p + -[H 1 (t)]- + -[H 1 {t)] + -p-+) 

+ ^r(p--p eq -(T)), 
J-i 



(35) 



(36) 



dp. 



dt 



-{E+-EJ)p-+ 



+ -(p__[/T 1 (t)]_+ -[#!(*)]__?_+ 

+ p-+[Hi(t)] ++ - [ff 1 (t)]_+p ++ ) - ip_ 

Li 



dp+_ dp*. 



dt 



dt 



(37) 
(38) 



where [H^t)]^ = {a\Hi(t)\(3) . Using Eq. Q, we can 
write the first two equations for the diagonal elements as: 



dp. 



dt 



= — ^ = ^(-r?A cosnt) (p_+ - p+_) 

+ - P- (P e + \(T) p e - 9 -(T))}, (39) 

J-i 



while the equation for the off-diagonal element p |_ is: 

fct = - S_)p_+ + ~(p__ - p ++ )(-77A cosOi) 



^P-+(-2r]Acosnt) - -^p_. 
n ±2 



(40) 



We look for a steady state solution of the form 

P—=r—, p++=r ++ (41) 

P-+(t)=r_+e im , p+_(t)=p* + (t), (42) 

where r aa / are complex constants. We find that in steady 
state the density matrix elements are given by the follow- 
ing expressions: 



1 

P— = o 



2 l+3(«,wo,r 2 ) x (7/J c ) 



1 

P++ = - 

p-+(*) - 



p? + (T)-l/2 



2 l + ff (n,w ,Ta)x(///c) 
ir 2 T?A 1 

2h [i + \T 2 (n -w )] 

p!?_(T)-^ + (T) 



(43) 
(44) 



l + 5 (ft,wo,T 2 ) x (7/J c ) 
P+-(*) = P-+(*)j 



exp(ifti) (45) 
(46) 



where V = p • £ ac /£, J// c = T 1 T 2 (t ? A /? 1 ) 2 /2, and 
#(fi, u; , T 2 ) = l/[l + (fi-u;o) 2 T|]. For a dipole moment 
p = 3.7 D, a large electric field £ ac = 3x 10 3 V/m, and 
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TLS energy splittings of the order of 10 GHz, the dimen- 
sionless factor r\ ps 0.005, and it decreases to a value of 
5 x 10 -8 when the amplitude of the applied electric field 
is £ae = 0.03 V/m. g(p,,wo,T2) is approximately equal 
to 1 when the ac driving frequency is resonant with the 
TLS energy splitting, i.e., SI ps uj . 

Notice that the off-diagonal elements of the density 
matrix are first order in r\ <C 1. They are oscillatory and 
small, as shown by the following numerical estimate. For 
large electric fields (£ QC = 3xl0 3 V/m), r) ps 0.005, T 2 = 8 
fjs at T = 0.1 A /h ps E/ft ps 10 GHz, and T x = 

8xl0~ 8 s, we obtain I/I c = T 1 T 2 (r?Ao/ft) 2 /2 ps 10 7 > 1, 
and pa 10~ 4 . For / < J C) we have \p-+\ ps 10~ 2 . 

Hence, the amplitude of the off-diagonal density matrix 
elements is very small for both />/ c and I <C I c - 

On the other hand, the diagonal elements recover their 
equilibrium values (pl 9 _(T) = exp(+E/2k B T)/Z and 
P++{T) = exp{-E/2k B T)/Z) in the limit of low electro- 
magnetic fields I <C I c . For large electric fields I ^> I c , 
they approach their steady state values p++ = p__ = 
1/2. As required, Tr{p} = p__ + p ++ = 1. In addition, 
the quantum expectation value of the z— component of 
the TLS spin is 



(S z ) = (p ++ - p--) 

tanh{E/{2k B T)) 



[l + (I/I e )xg(n,uo,T 2 )Y 



(47) 
(48) 



All the elements of the density matrix depend on the 
ratio I / I C (Q,T). We can write approximate expressions 
for (0__ and p ++ corresponding to the populations of the 
lower and upper TLS energy levels for both small and 
large electromagnetic fields. For I <C I c an d <7 ~ 1, we 
expand p and p++ to first order in I / I c to obtain: 



p__K P ^_(T){l-Lg)+\L g 



<■</ 



P+ + (T)(l--g) + --<?. 

J-C, 



1 I 
2~Z 



(49) 
(50) 



This 



means that unsaturated TLS have p++ ps and 

p ps 1. Thus the upper level is mostly unoccupied 

while the lower level is almost always occupied. 

On the other hand, for / I c , we can expand the 
steady state solution for p__ and p ++ given in Eqs. ([43]) 
and (gU to first order in (7// c ) _1 . The result is 

P-^l+ip^m-^ij-gy 1 (51) 

P++^~[p^-(T)-^](j-gy\ (52) 

Hence, for saturated TLS the lower and upper levels are 
almost equally populated, i.e., p++ ps p__ ss 1/2. This is 
to be expected since the TLS are constantly being excited 
by the ac electric field and de-excited by spontaneous 
and stimulated emission. Once the TLS are saturated, 
the populations of the lower and upper levels will have 
small deviations from their steady state values. We will 



look at the saturation effect in more detail in Section IVII 
where we plot the noise spectrum of a single fluctuator 
versus I / I c - This noise spectrum depends on the density 
matrix elements p ++ and p As one goes from the un- 
saturated regime to the saturated regime, the amplitude 
of p-- decreases by a factor of two. From Eq. ([24]) the 
polarization noise of a single TLS depends linearly on 

p Because p only decreases by a factor of 2 when 

the TLS are saturated, we will see that the saturation of 
TLS will not play an important role in the polarization 
and charge noise spectra. 



VI. RESULTS 

In this section we begin by studying the polarization 
noise spectrum of a single TLS fluctuator as a function of 
frequency and the electric field intensity ratio (I/I c ). We 
then obtain the total polarization noise by averaging over 
the distribution of independent two level systems. From 
this we get the charge noise that we will analyze at both 
low and high frequencies as a function of temperature 
and electric field intensity. 



A. Polarization Noise of One TLS Fluctuator 

Now that we have the matrix elements of the steady 
state density matrix in Eqs. (|4"3|) - (|46]) . we can use them 
to evaluate the expression for polarization noise found 
in Eq. (|24[) . In order to make further progress in eval- 
uating Eq. (|2"4|) , we need to know the polarization noise 
spectrum Sp (uj) of a fluctuating TLS. So we assume 
that a single TLS fluctuates randomly in time. Its elec- 
tric dipole moment fluctuates in orientation by making 
180° flips between 6 = 9 and (9 = 6 + 180°), resulting 
in a random telegraph signal (RTS) in the polarization 
Po(t) = ±pcos(#o)/V along the external field. Sp Q (u>) 
dominates at low frequencies when the system is in ther- 
mal equilibrium. So for Sp (u>), we will use a Lorentzian 
noise spectrum given by 32,33 



S Po (oj) = Awxw 2 



(1 



(53) 



where r is the characteristic relaxation time of the fluc- 
tuator, and w\ (W2) is the probability of being in the 
lower (upper) state of the TLS. Since the ratio of the 
probabilities of being in the upper versus lower state is 
wijw\ = exp(—E/ksT), and w± + u>2 = 1, the product 



AwiW2 



1 



cosh 2 (£/2fc B T) 



(54) 



So in Eq. (J24J), we replace Sp (cj) by the RTS noise spec- 
trum with a relaxation time t = T\ since this term is 
associated with a z , the longitudinal component of the 
TLS pseudospin. 
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FIG. 2: (Color online) Log-log plot of the low frequency po- 
larization noise of a single two level system given by Eq. (|56[) 
versus frequency for three different values of the TLS energy 
splitting as shown in the legend. For all 3 cases, T = 0.1 K, 
I/Ic = 1CT 4 , Q = E/h, and A = A = E/V2. 



FIG. 3: (Color online) Log-log plot of the high frequency po- 
larization noise of a single two level system given by Eq. (|56|) 
versus frequency for three different TLS energy splittings. 
The peaks appears when the resonance condition uj — E/h is 
satisfied. In all three cases, T = 0.1 K, I/I c = 10" 4 , Q = E/h 
and A = Ao = E/y2. No noticeable difference was obtained 
for I/I c = 10 7 > 1 (not shown). 



At high frequencies, Sp (u>±u>q) dominates and is asso- 
ciated with resonant processes. At high intensities when 
the driving frequency is close to the energy splitting 
ujo, saturation occurs, and the system is not in thermal 
equilibrium. So we will use 



S Po (lj±lj ) = 



(1 + (u±lj ) 2 t 2 ) 



(55) 



Since these terms are associated with the transverse com- 
ponent a x of the TLS pseudospin, we will set r equal to 
the transverse relaxation time T 2 . 

Putting this all together, we arrive at the following 
expression for the spectral density of polarization fluctu- 
ations of a single TLS: 



(^o 2 > 



E J l cosh 2 (E/2k B T) 1 
E 



(P~) 



1 + (u - E/h) 2 Ti 
T 2 



1 + (lj + E/hfTl 



dc term (56) 



The first term is the relaxation contribution and the mid- 
dle two terms are the resonance contribution. The dc 
term ensures that Sp,, (lj = 0) = 0. 

In Fig. [21 we show the spectral density Sp., (lj)/(SPq ) 
of polarization fluctuations of a single TLS at low fre- 
quencies. We consider three values of the TLS energy 
splitting as shown in the legend. In the low frequency 
range, the noise spectrum is dominated by the relaxation 
contribution which is a Lorentzian. Thus the noise spec- 
trum is flat for lj <C \/T\. As the frequency increases, it 
rolls over at lj ~ l/(Ti), and goes as 1/cj 2 for lj > \/T\. 
At very high frequencies (lj > \E/H\), it saturates to a 



constant (white noise) due to the resonance terms. In 
addition, we find that the low frequency noise is indepen- 
dent of electric field intensity ratio I / I c and the angular 
frequency fl of the ac driving field. 

The contribution of the resonance terms to the noise 
spectrum is negligible at low frequencies (lj -C E/h) 
as expected from simple numerical estimates. However 
these resonance terms become important at high frequen- 
cies as shown in Fig. [3] The peaks appear when the 
resonance condition lj = ujq = E/h is satisfied. The 
plot in Fig. [3] was obtained for the low value of the ra- 
tio I/I c = 10~ 4 . No noticeable differences were obtained 
when the ratio was increased to values as high as 10 7 (not 
shown) . 

Figured] shows the polarization noise power of a single 
two level system over a broad range of frequencies that 
covers both the resonance and relaxation contributions. 
At low frequencies there is a plateau. Between 100 kHz 
and 1 GHz, the noise spectrum decreases as 1/ui 2 due to 
the Lorentzian associated with the relaxation contribu- 
tion. At higher frequencies there is a resonance peak at 

lj = E/h. 

The effect of TLS saturation can be seen at high fre- 
quencies in the plot of the noise of a single TLS versus 
the ratio I/I c as shown in Fig. [5] We plotted the spectral 
density of polarization fluctuations of an individual TLS 
given by Eq. (|56p at a fixed high frequency (lj/2it = 9 
GHz) versus I/I c . Notice that the noise is constant as 
long as I <C I c , then decreases when the electric field 
intensity increases to a value comparable to the critical 
intensity, and reaches a value which is half the previous 
one for / 3> I c - This is in agreement with our estimates 

from Section fVl where we saw that as I / I c increases, p 

decreases by a factor of 2 from a value close to 1 to a 
value close to 0.5 when the TLS are saturated. 
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10 10 10 

Frequency [Hz] 

FIG. 4: Log-log plot of the polarization noise of a single two 
level system given by Eq. ([56} versus frequency for a broad 
range of frequencies. E = 5.2 GHz = 0.4 K, I/I c = 1CT 4 , 
Q = E/h, and A = A = E/y/2. 



obtain: 



(SPi) 



V 



A„ 



dA 



An 



dA P(A,A ) 



*Vr I 

E) L cosh 2 (E/2k B T)'l + io 2 T? 



Ao 
E 



(P— ) 



T 2 



1 + (w - E/K) 2 T 2 
T 2 



1 + (w + £yft) 2 T 2 



dc term 



V J(w; Q, T, J/J c ) 



2^2 

dc term, (57) 



where J(uj; f2, T, I/I c ) is the result of integrating over the 
distribution of TLS parameters A and Ao. 

To obtain the charge noise Sq(lu) from the polarization 
noise, we use Eq. (|27| with 




' 2p cos 6*c 

Vo 



(58) 
(59) 



We obtain 
Sq(w) _ 



4 , p 



(-y)T J(w; T, I II C ) + dc term. (60) 



3 y eL 



We can evaluate Eq. (fBUj) numerically using p = 3.7 D, 
P TiS = 10 45 (Jm 3 )" 1 , L = 400 nm, A = 40 x 800 nm 2 , 
V = AL, A mi „ = 0, A max /k B = A 0tmax /k B = 4 K, and 
A 0iro i„/fc B = 2x 10" 6 K. The results follow. We will 
drop the dc term from now on since it only affects the 
zero frequency noise. 



FIG. 5: Log-linear plot of the high frequency polarization 
noise for a single two level system given by Eq. (|56jl versus 
I/I c - The amplitude of the noise decreases by a factor of 
two as saturation is achieved for I/I c > 1. uj/2-k — 9 GHz, 
E/h = Q/2-k = 10 GHz, A = Ao = E/V2, and T = 0.1 K. 



Polarization and Charge Noise of an Ensemble 
of TLS Fluctuators 



Until now we have analyzed the contribution of a sin- 
gle fluctuator to the polarization noise spectrum. We 
now average the polarization noise of a single fluctuator 
over an ensemble of independent fluctuators in a volume 
V a . The distribution function over TLS parameters was 
given in Eq. ||5J) as P(A, A ) = Ptls/^-o- Using this, we 



Low Frequency Charge Noise 



Evaluating Eq. (|60p produces the normalized low fre- 
quency (uj -C \E/h\) charge noise spectrum Sq/c 2 shown 
in Fig. [HI It is flat at very low frequencies wT\ <C 1. As 
the frequency increases, it rolls over at uj w {Ti >max )~ l 
and decreases as 1// noise between approximately 10~ 3 
and 10 4 Hz. Here T\ max is the maximum value of T\. 
For T = 0.1 K, Ti lOT J ~ 10 4 s. 

We can obtain an approximate analytic expression for 
the low frequency noise from Eq. (|60[) in the follow- 
ing way. By low frequency, we mean that wr TO , n <C 1 
and ftw <C E. So we only keep the first term in 
J(w; Q, T, I/I c ) as defined by Eq. ([57]): 



5q(o;) 



dA 



dA P(A,A ) 



An 



L cosh 2 (P/2fc i3 T) J l + a; 2 T 1 2 



(61) 
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FIG. 6: (Color online) Log- log plot of the low frequency 
charge noise Sq/c 2 averaged over the TLS ensemble versus 
frequency for different temperatures. The symbols are the re- 
sults of numerically evaluating Eq. (|57[) with I/I c = 10 -4 and 
Q, = 10 GHz. The dashed lines are the low frequency approx- 
imation given by Eq. (|64[) . The solid line has a slope of —1 
corresponding to a perfect 1// spectrum. Between 10 -3 and 
10 4 Hz, the noise spectrum is very close to 1//. For example, 
in this frequency range for T = 0.1 K the noise spectrum goes 
as l/(/ 102 ). 



We change variables to the TLS energy splitting E and 
the relaxation time T\ . Using 



A 2 

E 2 



1 



Ii 



(62) 



we obtain 



Sq(lv) 



3W Vo 



1 - 



Ti 



dE P(E,Ti) 
Ti 



where P(E, T\) is given by Eq 
yields 

Sq(uj) ^ 4 , p s2 V Ptls 
e 2 ~ Z^eL' 2 



cosh 2 (E/2k B T) 1 + uj 2 T 2 
. Using (r TOi „/Ti) 



(63) 



< 1 



dE- 



deduced from current noise,™ The magnitude of this 
noise estimate is also in good agreement with our numeri- 
cal result from Eq. i.e., S Q (f = 1 Hz)/e 2 = 2 x 10~ 7 
Hz" 1 . 

We can also obtain this 1// noise result with the fol- 
lowing simple calculation. Consider an electric dipole 
moment p in a parallel plate capacitor at an angle #0 
with respect to the z— axis which is perpendicular to the 
electrodes that are a distance L apart. Assume that the 
electrodes are at the same voltage. When the dipole flips 
by 180°, the induced charges on the superconducting elec- 
trodes change from ^fp cos 6q/L to ±pcos9o/L. Let SQ 
denote the magnitude of the charge fluctuations. Then 
SQ = \2p cos 9 /L\. Hence the charge of the Josephson 
junction capacitor produces a simple two state random 
telegraph signal which switches with a transition rate 
t -1 given by the sum of the rates of transitions up and 
down 42, The charge noise spectral density is^ 3 - 



4» 



dih-tzy^-^iSQit^SQih)} 



(65) 



where the superscript i refers to the ith TLS in the sub- 
strate or tunnel barrier, and w\ (1U2) is the probability 
of being in the lower (upper) state of the TLS. In order 
to average over TLS, we recall from Eq. ^ that the dis- 
tribution function of TLS parameters can be written in 
terms of the energy and relaxation times) 2 ^ 



P(E,r) = 



Pt 



LS 



(2r V / l - Tmin/r) 



(66) 



At low frequencies ur mm <C 1, the main contribution to 
the spectral density comes from slowly relaxing TLS for 
which P(E, t) ~ Ptls/(2t). Therefore, the charge noise 
of an ensemble of independent fluctuators is 



dr 



(SQ 2 ) 



dE 



P 



TLS 



2t 



cosh I (E/2k B T) 1+u 2 t 2 



(67) 



n 



2tt. 



-V e 



co 2 T 2 

p \ 2 ksT 



Thus we obtain 1// noise that goes linearly with tem- 
perature. As Fig. [6] shows, this expression for the low 
frequency noise gives good agreement with our numerical 
evaluation of Eq. ([57)1 . Eq. ([M]) also agrees with Kogan 33 , 



cosh^ {E/2k B T) (SQ 2 ) is the square of the amplitude of charge fluctua- 
tions averaged over TLS dipole orientations. We assume 
that (SQ 2 ) is independent of E and r. At low tempera- 
tures (k B T -C E max ) and low frequencies (uT m i„ -C 1), 
(64) we recover Eq. (IM1) which describes 1 // charge noise that 
goes linearly with temperature. 

Still another way to obtain low frequency 1// noise is 
the following. At low frequencies the system is in thermal 
equilibrium, and we can use Eq. (|14[) from Section IIIII 



and Faoro and Ioffe^. To estimate the value of Sq/c 
from Eq. ([64]), we use p = 3.7 D, P TLS w 10 45 (Jm 3 )" 1 , 
L = 400 nm, A = 40x800 nm 2 , and V a = AL. At T = 0.1 
K and / = 1 Hz, we obtain S'g/e 2 = 2x 10~ 7 Hz -1 , which 
is comparable to the experimental value of 4 x 10~ 6 Hz -1 



Sq(uj) 2k B T tan S 
e 2 ~ e 2 /2C uj ' 



(68) 



The TLS contribution to the dielectric loss tangent tan 5 
was calculated by previous worker a 17 i 22 ' 23 who consid- 
ered fluctuating TLS with electric dipole moments. They 
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FIG. 7: (Color online) Log-log plot of the high frequency 
charge noise Sq/c 2 averaged over the TLS ensemble versus 
frequency. The contributions of the relaxation term and the 
resonance terms are plotted separately. The noise spectrum 
becomes white noise at frequencies ~ 10 MHz. The results 
were obtained using T = 0.1 K, O = 10 GHz, and I/I c = 10 7 . 




Temperature [K] 

FIG. 8: Plot of the high frequency charge noise SQ(u))/e 2 
averaged over the TLS ensemble versus temperature. Notice 
that the noise decreases gradually with increasing temper- 
ature. This plot was obtained for u) — 10 GHz, a driving 
frequency of Q = 10 GHz, and I/I c = 10 7 . 



found 



tan S — 



kP 2 Ptls 
6e' 



(69) 



where s' is the real part of the dielectric permittivity and 
Ptls the constant TLS density of states. By plugging 
Eq. KBSJ) into Eq. @SH and using Eq. |T3]) and V D = AL, 
we recover Eq. (|64[) . 



2. High Frequency Charge Noise 

Evaluating Eq. (|60[) numerically at high frequencies 
yields the normalized charge noise spectrum shown in 
Fig. [71 We have plotted the contributions coming from 
the relaxation term and the two resonance terms sepa- 
rately. As we described in the previous subsection, the 
relaxation term dominates at low frequencies and gives 
1// noise. In the high frequency range shown in Fig. [7j 
the relaxation term produces l// 2 noise. On the other 
hand, the resonance terms produce white (flat) noise that 
dominates at high frequencies. The two curves cross at 
approximately 10 MHz. 

Regarding the temperature dependence, we note that 
while the low frequency 1/ / noise is proportional to tem- 
perature, the high frequency white noise decreases grad- 
ually with increasing temperature as shown in Fig. [5] To 
understand this temperature dependence, note that at 
high frequencies the resonant terms dominate. These are 
the terms in Eq. (|57[) with (u) ± E/K) in the denomina- 
tor. The dominant contribution to the integral occurs 
at resonance (w = E/K) where the temperature depen- 
dence of the integrand goes as T 2 ~ 1/T and decreases 
as the temperature increases. However, this decrease is 
much stronger than seen in the ensemble averaged noise 
shown in Fig. [5] This may be because in obtaining the 
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FIG. 9: (Color online) Log-log plot of the charge noise 
SQ(ui)/e 2 averaged over the TLS ensemble versus frequency. 
Notice that the results for large and small values of I/I c over- 
lap. We have considered two temperatures, 0.5 K and 0.05, 
respectively, and the driving frequency was S7 = 10 GHz. 



ensemble averaged noise, one adds terms away from res- 
onance which have the opposite trend and increase with 
increasing temperature as T 2 _1 ~ T. 

In Fig. [S] we show our cumulative numerical results for 
the charge noise spectrum at both low and high frequen- 
cies. As we mentioned previously, there is no noticeable 
dependence on the electric field intensity ratio I/I c at 
either low or high frequencies. 

Our result of white noise at high frequencies dis- 
agrees with the experiments by Astafiev et al£^ who con- 
cluded that the noise increases linearly with frequency. 
However, we caution that the experiments were done 
under different conditions from the calculation. The 
experimentalists 21 applied dc pulses with Fourier com- 
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ponents up to a few GHz, possibly saturating TLS with 
energy splittings in this frequency range. They relied on 
a Landau-Zener transition to excite the qubit which had a 
much larger energy splitting, ranging up to 100 GHz, and 
measured the decay rate Ti of the qubit. Then they used 
Ti to deduce the charge noise spectrum at frequencies 
equal to the qubit splitting by using a formula 3 ^ derived 
assuming a stationary state. It is not clear whether it is 
valid to assume a stationary state in the presence of dc 
pulses which lasted for a time (~ 100 ps) comparable to 
the lifetime of the qubit. 

In contrast, in calculating noise spectra, we make the 
customary assumption of stationarity. We relate the 
charge noise spectra to the response to an ac drive for 
a broad range of frequencies. AC driving of qubits have 
been used in both theoretica l 35 ! 36 ' 37 and experimental 
studie d 38 ' 39 ' 40 ' 41 of qubits. It would be interesting to 
measure the frequency and temperature dependence of 
the charge noise in Josephson qubits in the presence of 
ac driving since no such measurements have been done. 



VII. SUMMARY 

To summarize, we have studied the effect of the satu- 
ration of TLS by electromagnetic waves on qubit charge 
noise. Using the standard theory of two level systems 
with a flat density of states, we find that the charge noise 
at low frequencies is 1// noise and is insensitive to the 
saturation of the two level systems. In addition the low 
frequency charge noise increases linearly with tempera- 
ture. As one approaches high frequencies, the charge 
noise plateaus to white noise with a very weak depen- 
dence on the driving frequency f2 and the ratio I/I c . 
We found that the high frequency charge noise decreases 
slightly with increasing temperature. 

Finally we note that while we have been considering a 
Josephson junction qubit, our results on charge and po- 
larization noise have not relied on the superconducting 
properties of the qubit. So our results are much more 
general and pertain to the charge noise produced by fluc- 
tuating TLS in a capacitor or substrate subject to a driv- 
ing ac electric field in steady state. 
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VIII. APPENDIX: DERIVATION OF THE 
CHARGE INDUCED ON THE ELECTRODES BY 
A DIPOLE 

In this appendix we derive Eq. (TTJ which gives the mag- 
nitude of the charge induced on the electrodes by a dipole 
Eq. (TTJ. Rather than use image charges which leads to 
an infinite sum, we follow Purcell4£ We start by con- 
sidering the simple problem of a charge q between two 



conducting plates connected by a wire so that they are 
at the same electric potential. The plates are parallel to 
the xy plane and separated by a distance L. The sum 
of the induced charges is —q. Let d u (di) be the perpen- 
dicular distance from the charge q to the upper (lower) 
plate. Let — q u (—qi) be the charge induced on the upper 
(lower) plate so that — q u — qt = —q. Notice that if the 
charge between the plates is doubled to 2q such that d u 
and di stay the same as before, the ratio q u /qi stays the 
same even though the total induced charge is — 2q. So let 
us replace q by a conducting plate with charge density a 
while still maintaining the distances d u and dg. The to- 
tal induced charge density is — 2a = —cr u — at where — a u 
and —a i are the induced charge densities on the upper 
and lower plates, respectively. Notice that 



<lu 



(70) 



Using Gauss' Law to find the electric field and the fact 
that the voltage difference between the middle plate and 
the upper plate equals the voltage difference between the 
middle plate and the lower plate yields 



at 



di_ 

d u 



(71) 



Now we return to the problem of the charge induced 
on the plates due a point charge q between the plates. 
From Eqs. (J7TJ) and ((701), we obtain 



q u = q- 



and 



d u 



(72) 



(73) 



Now suppose we have a dipole between the two con- 
ducting plates at the same potential. The dipole consists 
of two equal and opposite charges q+ = q and q- = —q 
separated by a distance d. The magnitude of the dipole 
moment is p = qd, and 9 is the angle between the z-axis 
and the dipole moment p. Let d u + {dz+) be the distance 
between q + = q and the upper (lower) plate, and let d u - 
{dg-) be the distance between g_ = — q and the upper 
(lower) plate. Then the charge induced on the lower plate 
is — qg where 



■ (d u+ - d u -) 
pcos 



(74) 



Similarly the charge induced on the upper plate is — q u 
where 



qu 



pcos 6 
L 



(75) 



Thus we recover Eq. (TTJ for Q, the magnitude of the 
charge induced on each electrode. 
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